What can we learn about cosmic structure from 

gravitational waves? 

Joan M. Centrella 

Laboratory for High Energy Astrophysics, NASA Goddard Space Flight Center, Greenbelt, MD 

20771 USA 


Abstract. Observations of low frequency gravitational waves by the space-based LISA mission 
will open a new observational window on the early universe and the emergence of structure. LISA 
will observe the dynamical coalescence of massive black hole binaries at high redshifts, giving an 
unprecedented look at the merger history of galaxies and the reionization epoch. LISA will also 
observe gravitational waves from the collapse of supermassive stars to form black holes, and will 
map the spacetime in the central regions of galaxy cusps at high precision. 


INTRODUCTION 

Most of the information we have about the universe in general, and cosmic structure 
in particular, has been gleaned through detections of electromagnetic radiation. For 
example, observations with optical telescopes have revealed galaxy morphologies and 
clustering patterns; studies of galaxy spectra uncovered the universal recession as well 
as the presence of dark matter. Radio telescopes have shown a variety of active galaxies 
with vast cosmic jets emanating from their centers. Observations of the microwave sky 
revealed the 3K background radiation and the imprints of structure at early epochs. 
And high energy detectors unveiled a variety of exotic phenomena, ranging from hot 
gas in galaxy clusters to accretion disks around massive black holes (MBHs) at the 
centers of galaxies. This wealth of this electromagnetic information comes generally 
from excitations in gases in the outer regions of stars, within galaxies and clusers of 
galaxies, and in accretion flows around compact objects such as black holes. 

As the 21 st century begins, a new observational window on the universe is be- 
ing opened by detectors designed to measure gravitational waves from astrophysical 
sources. Gravitational radiation is produced through the bulk motion of massive objects 
such as black holes, and reveals their dynamical behavior directly. This article begins 
with a short introduction to gravitational radiation, and then encompasses a brief tour 
of several gravitational wave sources that can contribute to our understanding of cosmic 
structure. 


A GRAVITATIONAL RADIATION PRIMER 

Gravitational waves are ripples in spacetime curvature that travel at the speed of light. 
They are typically generated by compact matter distributions that have time-changing 
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quadrupole moments' such as binary systems and nonspherical collapses, and carry 
information about the bulk motion of the sources. Since gravitational waves couple very 
weakly to matter, they easily escape from the centers of galaxies or collapsed objects, 
bringing information about these deep, hidden regions. Gravitational waves carry both 
energy and momentum. For a binary system, this means that both the binary separation 
and the orbital period will decrease due to the emission of gravitational radiation. 
Such behavior has been observed in the binary pulsar PSR B 191 3+16, providing strong 
confirmation of the existence of gravitational radiation as predicted by general relativity 
[ 2 ]. 

The characteristic dimensionless gravitational wave amplitude for a source of mass M 
located at a distance r can be estimated as (e.g., [3]) 
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Here, v ns is the characteristic nonspherical velocity in the source, Q is its (trace-free) 
quadrupole moment, and f? Schw its Schwarzschild radius. Thus, the strongest sources 
have large masses moving with v ns ~ c. For a stellar black hole system having mass 
M = 1OM 0 , this gives h ~ 10 -16 at distance r = 15kpc, and h ~ 10 -21 at r = 3000Mpc. 
The gravitational waves from a MBH system with M = 2.5 x lO 6 Af 0 at a distance 
r = 3000Mpc typically have amplitudes h ~ 10 -16 . 

The observed gravitational wave frequency of a source depends on its mass, physical 
scale, and redshift. For example, the nonspherical collapse of an object at redshift z will 
produce a burst of gravitational waves with a characteristic frequency 
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where t dyn = yj GMjR 3 is the dynamical timescale. The nonspherical collapse of a star 
at z ~ 0 with mass M ~ 1OM 0 near its Schwarzschild radius R ~ # Schw can produce a 
,burst of high frequency gravitational radiation at / burst ~ lkHZ. For the collapse of a 
supermassive star with M ~ 1O 6 M 0 at z ~ 2, the radiation is produced at low frequency, 
/burs. ~ 3.7mHz. 

A distorted black hole, such as might be formed from nonspherical collapse or the 
coalescence of two black holes or neutron stars, will emit gravitational waves as it 
settles down to a quiescent axisymmetric Kerr hole. These so-called ringdown waves 
have frequencies in the range [4, 5] 

/ring ~(1.2 — : 3.2)xlO^Hz(i^), (3) 

where the lower end of the frequency range corresponds to a nonrotating Schwarzschild 
hole, and the upper end to a maximally rotating Kerr hole. Thus, a distorted stellar 


1 Conservation of mass and momentum guarantee that there can be no monopolar or dipolar gravitational 

radiation, respectively. [1] 
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black hole with mass M ~ 2OA/ 0 in a nearby galaxy with z ~ 0 emits high frequency 
gravitational waves with / rjng ~ 1kHz. In contrast, a distorted MBH with mass M ~ 

2.5 x 10 6 M q rings down at low frequencies, with / ring ~ 4mHz at Z land / ring ImHz 
at z ~ 5. 

Another important class of gravitational wave sources consists of binary systems. The 
waves generated by a binary system whose components are on circular orbits and have 
comparable masses have frequency [6] 


, _ ^orb.tal TT 1 1 (™Y /2 
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where M is the total mass of the binary and a is the separation of the binary components. 
As the gravitational waves are emitted, the binary separation a decreases until it reaches 
a ~ 3 R Schw , by which point the components have begun to merge. This gives a frequency 
range for these inspiral waves [7] 

/inspiral ~ 400Hz ( (\+ z )M ) ' (5) 

Again we see that stellar black hole systems produce high frequency waves, so that 
/inspiral ~ 200Hz for a binary with total mass M ~ 2OM 0 at z ~ 0. MBH binaries 
emit waves at lower frequencies; a binary with total mass M ~ 2.5 x 1O 6 M 0 produces 
gravitational waves at frequencies / jnspira) < ImHz at z ~ 1 and /inspiral £ °' 3mHz at 
z ~ 5. 


DETECTING GRAVITATIONAL WAVES 

Broad-band gravitational wave detectors use laser interferometry to monitor the separa- 
tion between test masses located at a distance L from each other. When a gravitational 
wave impinges on a detector, the separation between the masses changes by A L, re- 
sulting in a strain amplitude h{t) = A L/L. Since gravitational wave detectors measure 
displacement h ~ 1/r (see Eq. (1)), improving the sensitivity by a factor of 2 increases 
the regions of the universe from which sources can be detected by a factor of 8. Elec- 
tromagnetic detectors generally measure energy flux ~ 1/r 2 , so a similar increase in the 
observable volume would require improving the detector sensitivity by a factor of 4. 

Ground-based kilometer-scale interferometers, such as LIGO, VIRGO, and GEO600, 
are sensitive to high frequency gravitational radiation in the range 10Hz < / Gw < 10 4 Hz. 
For example, the LIGO (Laser Interferometric Gravitational -wave Observatory) project 
[8] has 2 interferometers with arm length L = 4km, one in Livingston, LA and and 
the other in Hanford, WA, and a third with L — 2km at Hanford. LIGO has recently 
begun to take scientific data, and should be able to observe the dynamics of black holes 
with masses up to ~ few xlOM 0 as well as other high frequency sources. The first 
generation interferometers are expected to reach their best strain sensitivities h ~ 10 -21 
at frequencies near / ~ 200Hz. Second generation detectors can lower the broad-band 
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sensitivity to h ~ 8 x 10~ 23 and achieve even better sensitivity with “tunable” narrow 
band techniques [9]. 

The planned space-based LISA (Laser Interferometric Space Antenna) mission [ 10] is 
a joint NASA/ESA collaboration designed to detect low frequency gravitational waves 
in the range 10 _5 Hz < / GW < 1Hz. LISA consists of 3 spacecraft at the vertices of 
an equilateral triangle with arm lengths L ~ 5 x 10 6 km. The detector will follow 20° 
behind the Earth in its orbit about the Sun, and will observe gravitational waves from 
the dynamics of MBHs and other low frequency sources. LISA will reach its best strain 
sensitivity h ~ 10~ 23 around frequency / ~ 5mHz and should be capable of detecting 
MBH binaries to redshifts z > 10. 

The timing of millisecond pulsars provides a means of probing the very low frequency 
gravitational wave band, 10~ 9 Hz < / < 10~ 7 Hz. The rotation period of a millisecond 
pulsar is determined to high precision by measuring the arrival times of the pulses. 
Perturbations in the arrival times arise when the electromagnetic waves from the pulsar 
pass through distortions in spacetime caused by a stochastic background of gravitational 
waves [11], such as that from the ensemble of MBH binary coalescences across the 
universe. Current pulsar data sets reach a minimum characteristic strain h 4 x 10 -15 at 
frequencies / ~ few nHz; observations with a planned Pulsar Timing Array could lower 
this sensitivity to h ~ 10 -15 [12]. 


GRAVITATIONAL WAVES FROM MBH BINARIES 


MBH binaries are expected to result from the merger of galaxies [13, 14] and, in 
hierarchical structure formation scenarios, the merger of halos containing seed black 
holes [15, 16]. A variety of astrophysical mechanisms have been proposed to bring a 
MBH binary from an initally wide orbital separation a > lkpc to the regime in which 
gravitational radiation reaction will cause it to coalesce in less than a Hubble time. 
These include both dynamical friction from background stars [13, 17] and gas dynamical 
effects [18]. The time needed for a binary of total mass M and separation a to coalesce 
once gravitational radiation reaction becomes the dominant energy loss mechanism is 


[ 6 ] 
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assuming circular orbits and equal mass components. For example, a MBH binary of 
total mass M ~ 2 x 10 7 M Q starting from a separation a ~ 10 -2 pc will coalesce in a 
Hubble time ~ 10 10 years. 

The final gravitational radiation-driven coalescence can be thought of as proceeding 
in 3 stages: an adiabatic inspiral, followed by a dynamical merger, and a final ringdown. 
As the binary evolves from inspiral through merger and ringdown, the gravitational 
wave frequency increases. As long as any accompanying accretion disks are dynamically 
negligible, black hole binaries coalescing through the emission of gravitational radiation 
are solutions to the vacuum (i.e., source-free) equations of general relativity. As such, 
their timescales are proportional to the total mass and, after time scaling, all other 
properties of the dynamics and waveforms depend only on ratios involving the masses 
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and spins of the components. Thus, calcuations of the gravitational waveforms from all 
three stages of black hole binary coalescence can be easily scaled to apply to any masses 
and spins [3]. 

During the slow inspiral, the black holes are well -separated and spiral together adi- 
abatically. The resulting gravitational waveform is a “chirp,” which is a sinusoid that 
increases in frequency and amplitude as the orbital period shrinks; see Eq. (4). The 
gravitational waveforms from this inspiral phase are computed analytically using higher 
order post-Newtonian techniques in which the black holes are approximated as point 
masses [19]. These waveforms can be used as templates for data analysis algorithms 
based on matched filtering [20]. 

Typically, MBH binaries will be within the LISA frequency band for periods of 
roughly several months to years. If a sufficient number of cycles of the inspiral waveform 
are thus observed, both the so-called chirp mass J( c = (M,M 2 ) 3 / 5 /(M 1 + M 2 ) 1 / 5 and, 
with less accuracy, the reduced mass m {ed = A/,M 2 /(M, 4 -M 2 ) can be measured, as well 
as some information on the spins [21]. 

LISA’s ability to detect a binary black hole inspiral depends on both the mass and 
redshift of the system. In particular, the gravitational wave frequencies of a binary of 
masses M, and M 2 at redshift z are the same as those emitted by a binary local to the 
detector but having masses (1 +z)M l and (1 +z)M 2 [21]. Increasing the mass and/or the 
redshift, lowers the frequency; see Eqs. (3) and (5). For binaries with mass (1 + z)M ~ 
10 5 M q , most of the inspiral signal occurs in the frequency band in which LISA has very 
good sensitivity; this results in high signal-to-noise ratios (SNR) rsj 10 2 — 10 3 for such 
systems [21]. Overall, LISA should be able to detect the inspiral phase at SNR > 10 for 
equal mass binaries of total mass M ~ 1O 3 M 0 — 10 6 M 0 out to z ~ 9, and higher mass 
inspirals M ~ 1O 7 M 0 at lower redshifts z ~ 1 [21]. 

If MBH binaries proceed directly to the regime in which they coalesce under gravi- 
tational radiation reaction, recent calculations predict ~ 10 or more events per year will 
occur within LISA’s frequency band [22, 23, 24], Since many of these may originate at 
redshifts z>7, LISA can provide an outstanding probe of the high redshift universe, in- 
cluding the reionization epoch and the merger history of galaxies [24]. Observations in 
the nHz frequency range by pulsar timing arrays will detect the stochastic gravitational 
wave background from MBH binaries at lower redshifts, providing important constraints 
on that population [12]. 

The relatively slow inspiral is followed by the dynamical merger phase, in which the 
black holes plunge toward each other and merge into a single remnant, emitting a burst 
of gravitational waves. This stage is governed by strong nonlinear gravitational fields, 
and the resulting waveforms can only be calculated using full 3-D numerical relativity 
simulations [25, 26], Typically, LISA will be able to observe gravitational wave bursts 
from this dynamical merger phase on timescales of roughly minutes to hours, providing 
an outstanding probe of a regime that is expected to be phenomenologically rich. For 
example, the merger of two black holes may cause a substantial change in the direction 
of the spin axis of the more massive hole [27]. The jets emanating from the center of 
an active galaxy are believed to be directed along the spin axis of the central MBH 
[28]. Thus, if such a hole suffered a change in spin direction as a result of merging with 
another MBH, the change in direction of the jet would cause a new radio lobe to be 
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FIGURE 1. VLA radio observations of 4 radio galaxies showing the X-type morphology [30, 31, 32]. 
Reprinted with permission from Merritt, D. and Ekers, R. D., Science, 297, 1310 (2002). Copyright 2002 
American Association for the Advancement of Science. 


generated at some angle to the original one. This has been proposed as an explanation 
for the "winged" or "X-type" radio sources shown in Fig. 1 [29]. 

After the dynamical merger is complete and a common event horizon forms, the 
distorted remnant will ring down to a quiescent Kerr hole through the emission of 
gravitational radiation in quasi-normal modes. The ringdown waveforms are known 
analytically from calculations of black hole perturbations, and take the form of damped 
sinusoids [4]. LISA will observe these bursts on time scales of typically minutes to 
hours. Detections at SNR > 10 2 should be possible for masses M ~ 10 5 — 10 7 M Q out to 
z ~ 9 for any value of the spin [21]. For M ~ lO 8 A/ 0 at z ~ 1, the ringdown occurs at the 
edge of LISA’s low frequency sensitivity, and detection of the ringdown with large SNR 
requires large spin; c.f. Eq. (3) [33]. Observations of the gravitational waves from MBH 
ringdown will allow the identification of the mass and spin of the final Kerr black hole. 
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SUPERMASSIVE STARS AND MBH FORMATION 


The origin of MBHs is a topic of considerable interest for understanding the high redshift 
universe and the emergence of structure [34], Various scenarios have been proposed, 
including the collapse of baryonic gas clouds (c.f. [35]) and dark matter halos [36]. 
MBHs may form hierarchically through the coalescence of smaller “seed” black holes. 
They may also form through direct collapse of large gas clumps; in fact, a recent study 
reveals the intriguing possibility that such MBHs may form in binaries [37]. 

A large gas cloud undergoing direct collapse may first form a rotating supermassive 
star as an intermediate state. Such an object is dominated by radiation pressure and is 
subject to a general relativistic radial instability leading to further collapse [6]. Recent 
numerical simulations have demonstrated that a uniformly rotating supermassive star 
collapses coherently to a black hole containing. ~ 90% of the original mass; the rest of 
the gas forms a rotating disk outside the hole [38, 39]. A burst of gravitational waves 
from the collapse itself will be followed by ringdown radiation as the newly formed 
black hole settles into a quiescent state. For an interesting range of masses and redshifts, 
/ burst and / ring lie within LISA’s sensitive range; see Eqs. (2) and (3). If the supermassive 
star is differentially rotating, a non-axisymmetric instability could trigger bar formation 
before the catastrophic collapse; the rotation of this bar could generate quasi-peri odic 
gravitational waves in the LISA band even if a MBH does not form [40]. 


MAPPING MBH SPACETIMES 

The final class of gravitational waves sources that we consider here consists of compact 
remnants of evolved stars - white dwarfs, neutron stars, or black holes - orbiting close 
to the central MBH in a galaxy [41, 42], Compact remnants are expected to comprise 
a significant fraction of the total stellar population in the central cusp surrounding the 
MBH. If such an object is scattered into an orbit around the MBH, it will spiral in 
through the emission of gravitational radiation, eventually being captured by the MBH. 
If the remnant is on a tight eccentric trajectory, roughly the final year before capture will 
be spent executing ~ 10 5 — 10 6 orbits through the very strong field region near the event 
horizon of the MBH [43]. The gravitational radiation emitted during this stage can be 
detected by LISA at high accuracy [44, 45]. 

Since a stellar remnant of mass n spiralling into a MBH of mass M constitutes an 
extreme mass ratio binary, typically fx/M < 10~ 4 , the motion of the remnant can be 
treated perturbatively in the MBH spacetime. The emission of gravitational radiation 
as the particle moves in the MBH potential causes the orbit to evolve. Technical issues 
arise in computing the effects of radiation reaction on the orbits and the accompanying 
waveforms [45]. Recent work on nonequatorial circular inspirals around a massive Kerr 
black hole shows that the gravitational waveforms are modulated by the effects of MBH 
spin and frame dragging, as shown in Fig. 2 [43]. The modulation is stronger at later 
times when the remnant is closer to the MBH event horizon. Complicated behavior also 
arises from eccentric orbits [46] and the spin of the remnant. 

These trajectories and the resulting gravitational waveforms encode information about 
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Waveform at late times 



t x yu/M 2 


FIGURE 2. The gravitational waveform from an extreme mass ratio inspiral with /.i/M — 10 ~ 4 . The 
particle moves on a circular, nonequatorial trajectory around a central Kerr black hole which has nearly 
maximal spin, 0.998M. Reprinted from Ref. [43]. 


the multipolar structure of the gravitational potential of the central massive object. 
Measuring just three multipole moments can falsify whether this massive object is 
indeed a black hole [5]. LISA should be able to measure these extreme mass ratio 
waveforms out to a distance of a few Gpc if the remnant is a black hole with mass 
jU ~ 10 Mq, and out to a few x lOOMpc if the remnant is a white dwarf or neutron star 
[41, 42, 45]. Conservative estimates suggest that LISA will measure at least several such 
events per year, and possibly many more. 
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SUMMARY 


Gravitational wave detectors are poised to open a new observational window on the 
universe. Ground-based interferometers will observe the high frequency regime, 10Hz < 
/gw ~ 10 4 Hz, including sources such as the dynamics of stellar black holes. The low 
frequency region of the spectrum, 10Hz < / GW < 10 4 Hz, will be probed by the space- 
based LISA mission, and is especially important for understanding the emergence of 
structure in the universe. LISA will observe the coalescence of MBH binaries at high 
redshifts, yielding an unprecedented view of the merger history of galaxies and the 
reionization epoch. LISA will also be able to detect gravitational waves from the collapse 
of a supermassive star to a black hole, and will provide a high precision probe of 
the central regions of galaxy cusps. Pulsar timing arrays will observe the very low 
frequency regime, 10 -9 Hz < / < 10” 7 Hz, which includes the stochastic gravitational 
wave background from MBH binaries at lower redshifts. 
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